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Abstract 

We solved the three-nucleon (3N) Faddeev equation including relativistic features at incoming neutron 
lab energies El"-'' = 28, 65, 135 and 250 MeV. Those features are relativistic kinematics, boost effects and 
Wigner spin rotations. As dynamical input a relativistic nucleon-nucleon (NN) interaction exactly on-shell 
equivalent to the AV18 NN potential has been used. The effects of Wigner rotations for elastic scattering 
observables were found to be small. The boost effects are significant at higher energies. They diminish 
the transition matrix elements at higher energies and lead in spite of the increased relativistic phase-space 
factor as compared to the nonrelativistic one to rather small effects in the cross section, which are mostly 
restricted to the backward angles. 
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I. INTRODUCTION 

High precision nucleon-nucleon (NN) potentials such as AV18 CDBonn P], Nijm I, II and 
93 P] describe the NN data set up to about 350 MeV very well. When these forces are used 
to predict binding energies of three-nucleon (3N) systems they underestimate the experimental 
bindings of and ^He by about 0.5-1 MeV [J, |5|. This missing binding energy can be cured by 
introducing a three-nucleon force (3NF) into the nuclear Hamiltonian P|. 

The study of elastic nucleon-deuteron (Nd) scattering and nucleon induced deuteron breakup 
revealed a number of cases where the nonrelativistic description based on pairwise forces only is 
insufficient to explain the data. Generally, the studied discrepancies between a theory based on NN 
potentials only and experiment become larger with increasing energy of the 3N system. Adding 
now a 3NF to the pairwise interactions leads in some cases to a better description of the data. The 
elastic Nd angular distribution in the region of its minimum and at backward angles is the best 
studied example P, Q|- The clear discrepancy in this angular regions at energies below « 100 MeV 
nucleon lab energy between a theory based on NN potentials only and the cross section data can be 
removed by adding modern 3NFs to the nuclear Hamiltonian. Such a 3NF must be adjusted with 
each NN potential separately to the experimental binding of and ^He . At energies higher 

than ~ 100 MeV current 3NFs only partially improve the description of cross section data and the 
remaining discrepancies, which increase with energy, indicate the possibility of relativistic effects. 
The need for a relativistic description of 3N scattering was also raised when precise measurements of 
the total cross section for neutron-deuteron (nd) scattering P| were analyzed within the framework 
of nonrelativistic Faddeev calculations [lO||. NN forces alone were insufficient to describe the data 
above . .00 MeV, T,.e eSeCs d.e .o ■■e.a.ivi.ie .i„e,na.ie, considered in .ef. Q were oo^paraHe 
at higher energies to the effects due to 3NFs. These indications show the importance of a study 
taking relativistic effects in the 3N continuum into account. 

The estimation of relativistic effects on the binding energy of three nucleons has been the focus 
of a lot of work. Basically two different aproaches have been followed: one is a manifestly covariant 
scheme linked to a field theoretical approach 111 , Il3, Im , the other one is based on relativistic 
,„aotu. mechanics for.n.ted on sp.el.e .^e.nrW . M.n.owsH .pace Q Q Q Q- 
Within the second scheme the relativistic Hamiltonian for on the mass shell particles consists of 
relativistic kinetic energies and two- and many-body interactions including their boost corrections, 



which are dictated by the Poincare algebra 



Ly-body mterac 
P, 15, 16,13 



The applications of these two types 



of approaches to the 3N bound state have led to contradictory results. In the approach based on 
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field theory, relativistic effects increase tlie triton binding energy, while in the approach based on 
relativistic Hamiltonians they decrease the triton binding energy. This requires further insights 
which is beyond the scope of this paper. 

Due to the increased complexity of 3N scattering calculations as compared to the bound state 
problem no results for the 3N continuum including relativity are available. In order to extend the 
Hamiltonian scheme in equal time formulation to 3N scattering one needs as a starting point the 
Lorentz boosted NN potential which generates the NN t matrix in a moving frame via a standard 
Lippmann-Schwinger equation. Such potentials have been worked out and applied to a 3N bound 
state in 1131. The results obtained supported the relativistic effects found before in a relativistic 

n 

quantum mechanics approach [19|]. The starting point for a NN potential in an arbitrary moving 
frame is the interaction in the two-nucleon cm. system, which enters a relativistic NN Schrodinger 
or Lippmann-Schwinger equation. It differs from the nonrelativistic Schrodinger equation just by 
the relativistic form for the kinetic energy. The current realistic NN potentials are defined and 
fitted in the context of the nonrelativistic Schrodinger equation. Up to now NN potentials refitted 
with the same accuracy in the framework of the relativistic NN Schrodinger equation do not exist. 
In such refitting was omitted and an analytical scale transformation of momenta which relates 
NN potentials in the nonrelativistic and relativistic Schrodinger equations in such a way, that 
exactly the same NN phase shifts are obtained by both equations, was employed j2o|. 

Though this transformation is not a substitute for a NN potential with proper relativistic 
features it can serve as a first ste p to illustrate the effects of Lorentz boosts on NN potentials. 
Such an approach was applied in 1^ and we also will follow it in the present study to get the 
first estimation of relativistic effects in the 3N continuum. In this first study we would like to find 
out what are the changes of elastic nd scattering observables when the nonrelativistic form of the 
kinetic energy is replaced by the relativistic one and a proper treatment of boost effects and effects 
due to Wigner rotations of spin states is performed. 

The paper is organized as follows. In Sec. II we lay out the relativistic features underlying our 
treatment for a relativized Faddeev equation in the 3N continuum. This incorporates the definition 
of the boosted two-body force, the various two-and three-body states in general frames, the Wigner 
rotations and the singularity structure of the relativistic free 3N propagator. Our manner to treat 
the 3N Faddeev equation is guided by the lines presented in [2]] for the nonrelativistic case. In Sec. 
Ill we focus on the relativistic NN potential and discuss the quality of different approximations for 
the boosted potential. As a consequence, in this first study we restrict ourselves to the leading order 
relativistic term in the expansion of the boosted potential. In Sec. IV we apply our formulation 
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based on a relativistic NN interaction which is exactly on-sheh equivalent to the nonrelativistic 
AV18 potential and solve the relativized 3N Faddeev equation with different approximations for 
the boost. We show and discuss results for elastic Nd scattering. Sec. V contains a summary and 
outlook. 



II. FORMULATION 



The nucleon-deuteron scattering with neutron and protons interacting through a NN poten- 
tial V alone is described in terms of a breakup operator T satisfying the Faddeev-type integral 
equation |21 



T\^ > = tP\^ > +tPGoT\^ > . (1) 

The two-nucleon (2N) t-matrix t results from the interaction V through the Lippmann-Schwinger 
equation. The permutation operator P = P12P23 + P13P23 is given in terms of the transposition 
Pij which interchanges nucleons i and j. The incoming state \(f) >= \qo > {(pd > describes the 
free nucleon-deuteron motion with relative momentum qq and the deuteron wave function \(j)d >■ 
Finally Go is the free 3N propagator. The physical picture underlying Eq.(P) is revealed after 
iteration which leads to a multiple scattering series for T. 



The elastic nd scattering transition operator U is given in terms of T by 
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U = PGq^ + PT. (2) 

This is our standard nonrelativistic formulation, which is equivalent to the nonrelativistic 3N 
Schrodinger equation plus boundary conditions. The formal structure of these equations in the 
relativistic case remains the same but the ingredients change. As explained in the relativistic 
3N Hamiltonian has the same form as the nonrelativistic one, only the momentum dependence of 
the kinetic energy changes and the relation of the pair interactions to the ones in their corresponding 
cm. frames changes, too. Consequently all the formal steps leading to Eqs.(^ and @ remain the 
same. 

The relativistic kinetic energy of three equal mass nucleons in their cm. system can conveniently 
be presented by introducing the free two-body mass operator. Let k and —A; be the momenta in 



one of the two-body subsystems, then 2uj{k) = 2y m? + k is the momentum dependent 2N mass 
operator and the 3N kinetic energy can be written as 



Ho = M2uj{k)Y + q^ + Jm? + q^, (3) 



where q is the momentum of the third particle and —q the total momentum of the chosen two-body 
subsystem (m is the nucleon mass). Any two-body subsystem can be chosen. As introduced in 
[2^ the pair forces in the relativistic 3N Hamiltonian living in moving frames are chosen as 

V{q) = ^{2uo(k)+vf + q^ - ^ {2Lo{k)f + q^ , (4) 

where V{q) for g = reduces to the potential v defined in the 2N cm. system. Note that also 
in that system the relativistic kinetic energy of the two nucleons has to be chosen, which together 
with V defines the interacting two- nucleon mass operator occuring in Eq.(|H). 

Let us now firstly regard the 2N subsystem. The standard nonrelativistic 2N Lippmann 
Schwinger equation turns now into a relativistic one, which in a general frame reads 

t{k,k'-,q) = V{k,k';q)+ [ d'k" ^ V{k,k " ■,q)t{k " ,k '-q) 



We refer to k ']q) as the boosted 2N t-matrix like we talk of the boosted 2N potential in 
Eq.dll). 

Using @ the relativistic 2N Schrodinger equation for the deuteron in a moving frame can be 
cast into the form 

<t>d{k) = , \ ^ / d^'k'Vik^k'-mdik'), (6) 



^Mj + q,' - ^{2oo{k)f 



where yMj + q^' is the energy of the deuteron in motion and M^t its rest mass. This equation is a 
good check for the correct numerical implementation of the boosted potential V{q) as will be used 
below. 

The new relativistic ingredients in Eqs.(^ and Q will therefore be the boosted t-operator and 
the relativistic 3N propagator 

where Hq is given in Eq. (3) and E is the total 3N cm. energy expressed in terms of the initial 
neutron momentum ^0 relative to the deuteron 



E = ^{M,f+q^^ + ^m^ + q^\ (8) 
Currently Eq.(^ in its nonrelativistic form is numerically solved for any NN interaction using a 



momentum space partial wave decomposition. Details are presented in ref. 



2l|. This turns Eq.© 



into a coupled set of two- dimensional integral equations. As we show now, in the relativistic case 



we can keep the same formal structure, though the kinematics and the momentum representation 
of the permutation operator P is more complex and boosted t-operators as well as Wigner rotations 
will appear. 

In the nonrelativistic case the partial wave projected momentum space basis is 

\pq{ls)j{\^)IJ{t^)T >, (9) 

where p and q are the magnitudes of standard Jacobi momenta (see and {ls)j two-body 

quantum numbers with obvious meaning, (A1/2)I refer to the third nucleon (described by the 
momentum q), J is the total 3N angular momentum and the rest are isospin quantum numbers. 
This is now to be generalized to the relativistic case. 

We regard firstly the two-nucleon system and replace the nonrelativistic relative two-nucleon 
momentum p by k, where k and —k are related to general momenta of two nucleons, say p2 and 
P3, by a Lorentz boost: 

k = k{p2,P3) 

= ^{P2-P3-P23 , ^ ^^ ), (10) 

E23 + V -E'23 - 'P23 



with Ei = ^Jrrfl + Pj^, E23 = E2 + E^ and P23 = P2+P3- This is a relativistic generalization of the 
nonrelativistic relative momentum p. 

The individual momentum state with momentum k for an on-the-mass-shell particle with mass 



m is defined in terms of the state in the rest frame as 
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25| 



u{m)\o^^> = ^^^\kfi>. (11) 

Here U{P{k)) is the unitary operator related to the special (along k) boost matrix (3{k), ^ the spin 



magnetic quantum number and k = {ujm{k), k) with uJm{k) = \J m? + k"^. Note that by definition 
II does not change. The 4x4 matrix (3{k) is given in the Appendix. 

A following unitary operation related to the general boost A leads to the well known Wigner 



rotation 
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U(K)\k^i> = /^^[/(A)C/(/3(A;))|0/z > 

V ^m{k) 



= U{mmf3-\k')Ap{k))\6f^ > (12) 

V ^mik) 

The argument in the second unitary operator is a rotation matrix 

RiA,k) = (3-^{k')A(3{k) (13) 



and often denoted as Wigner rotation. Further the on sheU momenta k and k' are related by 
k' = Ak. 

i?(A, k) being a rotation yields 



U{A)\kfi > = um'))Y,D}^{R{A,k)W > J^^- (14) 



and using (|TT|) again to 

U{A)\kfi > 



^^J2Dl,^{R{A,km'^^'>. (15) 



Here D^,^ are the standard SU(2) Wigner D-matrices 1.26'] and their arguments Euler angles which 
are related to R as shown below. 

Now we apply a general boost to the noninteracting two nucleon state 

|^;0/X2/i3> = \kfJ,2,-kn3 > ■ (16) 

In the notation to the left we changed from the two individual momenta to the relative momen- 
tum k and the total two-nucleon momentum zero. Since the two nucleons are noninteracting U (A) 
acting on a two-body system is a tensor product 

UiA) = U2iA)U3{A) (17) 

acting on the two spaces of nucleons 2 and 3. We choose A as P{P), where P = {Pq,P) = 
{ujm{p2) + Wm(p3),P2 + Ps)- That boost matrix (3{P) maps k into p2 and —k into p^. We obtain 

C/(/3(P))|fc;0/i2/i3 >= U2{l3(P))\k^,2 > U3{/3{P))\ - ^/xs > 



/u^^^l ^^(^^^^^^^))|.^, > J^E4m3(^(^(^)'-^))I^^3^;, > . (18) 

V ^m[K) ^/^ V ^m{K) ^i^ 

Of course p2 and p^ are given by the inverse relation to (10) and the related expression for —k. 

On the other hand the noninteracting two-nucleon system with total momentum zero can be 
considered as one object with a mass Mq = 2ujrn{k) and therefore according to the general relation 
((TT|) one obtains 



i7(/3(P))|fc;W3> = y^^^^5^|fc;P/i2/X3>, (19) 



and we end up with 

d{p2 P3 
d{P k) 



k-Pf,2f^3> = I ^^74^ |5^Z);,^^(i?(/3(P),A;)) 



^2 



where 



Y,Dl,^JR{PiP),-k))\p2^^',ps^i's >, (20) 
M3 



d{p2 P2) I Mo UJmim) UJmim) _ ,.2/- 



4^^N\p2,P3) (21) 



5(P A;) WAfo(^) Wm(A;) ^^(fc) 

is the Jacobian for the Lorentz transformation from (pg, Ps) to (P,k). 

This relation generaHzes the one used in jl^] and [jj] from the spinless case to the one with 
spin. In the nonrelativistic case the D-matrices reduce to Kronecker symbols and the Jacobian 
is one. Thus |p2P3/^2/^3 > equals directly |/cP/i2/U3 >, where k equals the nonrelativistic relative 
momentum. 

The next step is the transition to partial waves. Firstly one defines the two-body orbital angular 
momentum states 

\k l;0fiifi2f^3 > = / dkYl^,{k)\k;6fi2f^s > . (22) 

Then we couple with the total spin s in the two-nucleon cm. system to the total angular 
momentum j and its magnetic quantum number /x: 

^11 

\{ls)jkfi;0> = ^ ^(-/i2-/i3|s/Xs)(///iS^s|j»|Hw/^2At3;0 > . (23) 

The special boost (3{P) leads then to 



U{m)mjkfi;6> = ^'^^^^\{ls)jkf,;P> (24) 



and applied individually leads finally to 



\ils)Jk^,■p> = I 1^ Yl E(^/^4''3k/i.)(/w5/x,|j^) / dkYl^^ik) 

X E^iM2(^(/^(^)'^))E^4M3(^(^(^)'-^))l^"*2/^2P3^3 > • (25) 
^2 Ms 

This is the connection of the partial wave projected two-nucleon state with internal momentum k 
and total momentum P to arbitrary individual momentum and spin states. 

Another requisite is the determination of the Euler angles (a,/3, 7) in the spin 1/2 D-matrices. 
According to (13) the two 4x4 matrices 

RiPiP),k) = (3-\p2)(3{P)f3{LOm{k),k) 
RiP{P),-k) = (3-\p3)PiP)PMk),-k) (26) 



representing rotations have the structure 

VO M J 

where M is the unitary 3x3 matrix for the Wigner rotation. It has generahy the form 
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/ cos a cos f3 cos 7 — sin a sin 7 sin a cos f3 cos 7 + cos a sin 7 — sin f3 cos 7 \ 
M = — cos a cos (3 sin 7 — sin a cos 7 — sin a cos /? sin 7 + cos a cos 7 sin /? sin 7 . (27) 
\ cos Q sin P sin a sin /3 cos f3 / 

This determines the 3 Euler angles. The matrix M related to the first equation in (|26j) is given in 
the Appendix. 

It remains to add the third free particle whose momentum pi together with the total two-nucleon 
momentum P adds up to zero in the 3N cm. system. 

Sticking to our standard nonrelativistic notation we denote the orbital angular momentum of 
that third particle by A and couple it with its spin to its total angular momentum /. Then the 3N 
partial wave state is 

\kq = pia > = \kpi{ls)j{X^)I{jI)JM > \{t^)TMT > 

11 1 

= N{p2,P3) i2^^2^^J'3\s^^s){l^^lS^^s\j^J'){)^fJ'\^^J'l\I|J'I){j|J'I^J'I\JM) 

At2/i3Ms |lll^'2^l's MlA^AM/M 

dpiY,\{pi) I dkYlik)Dl(R{l3{P)M Dl(R[f3{Pl-k)) 



|P2At'2#3/^'3Pl^l > • (28) 

In that expression p2 and are functions of k and q = —pi . 

For the evaluation of the partial wave representation of the permutation operator P we need 
the projection of that state \kq = pia > onto < pi ' ^'iP2 'fJ'2P3 Val- Doing that one encounters 
(P = -Pi) 

S{P2' -P2{k,P))m' -P3ik,P)) 

= 7^^T^^s{k - k{p2 ',m 'MP -n'- P3 ') (29) 
I d{k, P) I 

This is verified for instance by integrating both sides over p2 ' and ' and by converting the 
integral on the right hand side to an integral over k and P. Thus using (28) one obtains 

1 



< pimip2m2P3mz\k q a > = S{pi+p2+p3)- ^ 

N{P2,P3) 

S{q -Pi) ^{k - k{ p2,P3 )) 
q pi k k{ P2,P3) ) 
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(ll^2^l^3\sf^s){lfJ'iSfis\jfJ')i>^f^x^rni\IiJ,i){ji^Il^i\JM)Yi^^{p^ )y^^(A:( P2,P3) ) 



~2 2 



Dl,^,iRiPiPi P2,P3) ),k{ P2,P3) )) 

Di,^,{R{PiPi P2,P3) ), -k{ P2,P3) ))■ (30) 

This is the basic expression needed for the evaluation of the partial wave representation of 
the permutation operator P. Equipped with that, projecting Eq.(l) onto the basis states (28) one 
encounters like in the nonrelativistic notation [2J| 

1 < kqa\P\k'q'a' >i = i < kqa\k'q'a' >2 +i < kqa\k'q'a >3= 2 i < kqa\k'q'a >2 • (31) 

This is evaluated by inserting the complete basis of states \pifiiP2fJ'2P3fJ'3 > and using (30). The 
result is worked out in the Appendix. It can be expressed in a form which resembles closely the 

"1,3 



one appearing in the nonrelativistic regime 



^<kqa\P\k'q'a'>^ = //-^^^ " 

1 1 



-1 A;'+2 A;'''+2 

1 1 

Gaa'{q,q',x), (32) 



Niiq,q',x) N2iq,q',x 

where all ingredients are given in the Appendix. 

It remains to regard the free propagator adjacent to the permutation operator P in Eq.(^. 
Since only momenta are involved in the propagator the convenient formal steps outlaid in 27] can 
be shown in a momentum vector notation thereby simplifying the notation. Let \k,q>i denote the 
3N state expressed in vector momenta analogous to (20) and neglecting spin and isospin degrees 
of freedom. The index 1 indicates as above that the 2N subsystem (23) has been chosen which 
is described by the internal momentum k. Similarily an index 2 indicates the choice of the (31) 
subsystem. Then one obtains 

1 < kq \G0P12P23\k 'q' >i=i< kq\GQ\k 'q' >2 
= j Wdpi 1 < kq\Go\pip2P3 > ^iPi +P2+P3) < PiP2P3\k 'q' >2 

= I Wdpi . ^, , 1 < kq\pip2P3 > 

•' E + ie — J m? + p^ — J m? + P2 — J w? + p^^ 

5{pi +P2+P3) < PiP2P3\k 'q' >2 

\{dpi- ^ / , ■^==5{q-pi)5{k-k{p2,P3)) 

E + le — Jm^ + p^ — J + P2 — ^Jm^ + P3 

^(^(pi +P2 +^3)^(9 ' -P2)^(A: ' - A;(p3,Pi))- 



N{p2,P3) N{ps,pi) 



E + ie — \J + g ^ — \J m? + q '"^ — yjrn^ + (^+ 9*') 



==^S{k - k{q',-q- q")) 
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— —rS(k'-k(-q-q',q])-—^^—— (33) 
N{q',-q-q') ^ ^ ^ ^ ^^'> N{-q - q ' ,q) ^ > 

We recognize that the free propagator depends on q, q' and x = q - q Hke in the nonrelativistic 
case. After partial wave decomposition there arises an integration over the interval [—1,1] for 
X. This leads to the well known logarithmic singularities, which have been well studied for the 
nonrelativistic free propagator. That nonrelativistic propagator results simply by expanding the 
square roots in Eq. (|33)) and keeping the leading terms. Like in the nonrelativistic case it is now 
convenient to put the free propagator into the form oc l/(xo — x + ie) . A simple algebra leads in 
obvious notation to 

1 A 



E + ie - E^- > - -Eg+g ' xo - X ' 

with 



(34) 



A 



E — E^ — Eg' ' + Eg-^g- 



2qq' 

'2qq' 

Altogether we end up with the infinite system of coupled integral equations analogous to the 



{E-Eg^-Eg^,r-m'-q^-q'^ 
^0 = ^-1^. • (35) 



one in the nonrelativistic case 



21 



221: 



< kqa\T{E)\(j) > = < kqa\tP\(j) > +yy^ / dq'q'^ / dx 



X 



G^^>{q,q',x) <n2q'a'\T{E)\^> 



Ni{q,q',x) N2{q,q',x) ttz'-' 

X . (36) 

xq + ie — X 

The geometrical coefficients Gaa'iQ,Q' the coefficients Ni{q,q',x) and N2{q,q' ,x), and the 
momenta vri and tt2 stem from the matrix element < kqa\P\k'q'a' > of the permutation operator 
fEas. UCGf) and (jCTf) ). The quantum numbers in the set a differ from those in a only in the orbital 
angular momentum / of the pair. 

As mentioned earlier the main problem of treating Eq. H36|) is caused by the singularities of the 
free propagator Gq which occur in the region of q and q' values for which |xo| < 1. In addition 
at g = go there is the singularity of the 2N t-matrix which occurs in the ^Si Di partial wave 
state, where the deuteron bound state exists. The method to treat this singularity is described in 
detail in ^|. It amounts to separate all channels a which are "deuteron" -like (angular momentum 
quantum numbers / = or 2, s = 1 and j = 1) from the others. In the "deuteron" -like channels one 
separates the bound state pole and treats it by subtraction Similarly, the treatment of the Gq 
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singularities follows the nonrelativistic case as described in detail in ref. |2]|. The only difference 
is that the relativistic boundaries of the q,q'-values at which |xo| = 1 differ from the nonrelativistic 
boundaries. Equation ()36p is solved by generating its Neumann series, which is then summed up 
by the Pade method. 

Due to short-range nature of the NN force it can be considered negligible beyond a certain value 
jmax of the total angular momentum in the two nucleon subsystem. Generally with increasing 
energy jmax will also increase. For j > jmax we put the t-matrix to be zero, which yields a finite 
number of coupled channels for each total angular momentum J and total parity vr = ( — j'^'^ of 
the 3N system. To achieve converged results at our energies we used all partial wave states with 
total angular momenta of the 2N subsystem up to jmax = 5 and took into account all total angular 
momenta of the 3N system up to J = 25/2. This leads to a system of up to 143 coupled integral 
equations in two continuous variables for a given J and parity. 



III. THE BOOSTED POTENTIAL 



As dynamical input we used a relativistic interaction v, which is defined as partner of the 
relativistic kinetic energy, generated from the nonrelativistic NN potential AV18 according to the 
analytical prescription of ref. [3]. For the convenience of the reader we repeat the main points 
of this transformation. Having a NN potential f "'' which provides a nonrelativistic t-matrix t"'^ 
obeying the Lippmann-Schwinger equation with the nonrelativistic form of the free propagator 
one can apply an analytical transformation of momenta to obtain an exactly on-shell equivalent 
relativistic potential f'"' which provides the corresponding relativistic t-matrix f'. This t-matrix 
obeys the Lippmann-Schwinger equation with a relativistic form of the free propagator. This 
analytical transformation for the potentials is [2^ 

^ ^ 1 — — 1 

''^^'^'^ = ^(m''^"''^^^"^'^"''^/^ ^^^^ 

where 



knr =1 knr \ = v2m]j \J k + m? — m 

h{knr) 



h2 / h2 



In ref. |l^ it was shown that the explicit calculation of the matrix elements V{k,k';q) according 
to Eq. (4) for the boosted potential requires the knowledge of the NN bound state wave function 
and the half-shell NN t-matrices in the 2N cm. system. In this first study we do not treat the 
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boosted potential matrix element in all its complexity as given in ref. but restrict ourselves to 
the leading order term in a qjuj and vjio expansion 

Vi^S'-A) = v{k,k') 

X [ 1 - J' ^= ]. (39) 

S\Jm? + k '^\jm^ + {k'Y 

It is therefore important to check the quality of such an approximation. To that aim we 
calculated the deuteron wave function <j)dik) for the deuteron moving with momentum q using 
Eq.@. This wave function depends only on the 2N cm. relative momentum k inside the deuteron 
and is thus independent from the total momentum q. 

We show in Fig. ^ the binding energy and the D-state probability Pd defined through 
Eq. @ as a function of the initial nucleon lab. energy using the approximation given in Ea. H39() . 
In addition, the results for two more drastic approximations are given. In the first one the boost 
effects are neglected completely 

V{k,k';q) = v{k,k'), (40) 

and in the second one the k-dependence of the first order relativistic correction term is omitted 

V{k,k';q) = v{k,k') (^1--S^y (41) 

When the boost effects are fully taken into account the solution of Eq. © must provide exactly 
the deuteron binding energy and the D-state probability equal to the values for the deuteron at rest. 
It is seen in Fig^that neglecting the boost totally or omitting the k-dependence of the first order 
term is a poor approximation, especially at the higher energies. In contrast, the approximation 
given in Eq.ljSH) appears acceptable, even for the strongest boosts. Relying on that result we have 
chosen the expression p9|) for the boosted potential in the following investigations. 



IV. RESULTS 



The solution of the 3N relativised Faddeev equation including Wigner spin rotations increase the 
computer time drastically. This is caused by the calculation of the permutation matrix elements 
for the high partial waves. Therefore to study the effects of the Wigner spin rotations on the 
elastic scattering observables we restricted ourselves to the j < 2 partial wave states. We checked 
that to get converged results for the permutation matrix elements one has to take into account the 
expansion coefficients a^^^itj^f^^^^ {q, q') of Ea. ()C4|) with L,L' up to L,L'< 2. Those coefficients 
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FIG. 1: The deuteron binding energy Ed and the D-state probability Pjj. The values for the deuteron at rest 
are given by the solid horizontal lines. At different incoming nucleon lab. energies, related to the relative 
momentum qq in Eq. (6), the approximations given in Eas. H40(l . (|41(l . and H39II provide results which are 
shown as the dashed-, dotted-, and dashed-dotted lines, respectively. The calculations have been done with 
the AV18 potential. 

were obtained by numerical integrations over the directions q and q' with 23 gaussian points for 
the polar and azimuthal angles. We found that the changes of the cross sections due to Wigner 
spin rotations are small and stay under 1%. For spin observables these changes are slightly larger 
but they do not exceed 5% with the exception of angular regions around zero crossings and small 
values of the observables. Thus when performing the fully converged calculations with j < 5 and 
J < 25/2 we neglected the Wigner spin rotations completely. This might be different for breakup 
observables, which deserves another investigation. 

In FigOwe show our results for the nd elastic scattering cross sections at four energies together 
with experimental pd data. In addition to the nonrelativistic prediction, based on the solution of the 
3N Faddeev equation with the nonrelativistic form of the free propagator Gq and partial wave states 
constructed with standard Jacobi momenta, also our relativistic results using the approximations 
according to the EQs. (|in|) . (|1T|). and for the boosted potential matrix elements are presented. 
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FIG. 2: The differential cross sections for Nd elastic scattering at various energies. The soUd hne is the 
result of the nonrelativistic Faddeev calculation with the AV18 potential. The relativistic predictions based 
on the approximations H4U|) . H41|l . and H39|) for the boosted potential keeping all other relativistic features 
unchanged are shown by the dashed, dotted, and dashed-dotted lines. The pd data at 28, 65, 135 and 
250 MeV are from ref. 

0, 0> 0. and 0, 

respectively. 

Thereby the other relativistic features in Eq. (36) have been kept. It is only the total neglection of 
the boost effect in V{q) which leads to a clearly visible deviation from the nonrelativistic results. 
Taking the boost effect into account according to the approximations (41) and (39) reduces the 
effect drastically and only at the largest angles deviations from the nonrelativistic results are 
discernible. This is better seen in Fig. 3, where the quantity 



/ da_\rel ( da \nrel 



' da \nrel 

^dn) 



(42) 



expressed in percentage is shown for the three approximations. Thus significant effects of relativity 
occur at higher energies and they are restricted to the backward angle region (^c.m ^ 160°). They 
increase the nonrelativistic cross sections by up to 2%, 6%, 5 — 15%, and 10 — 23%, for 28, 65, 
135, and 250 MeV, respectively. For 6c.m < 160° the effects of relativity are much smaller. At 250 
MeV where they are largest, they increase the nonrelativistic cross section by no more than ~ 5% 
in the minimum around 9c.m = 130°. At forward angles the largest effects are at Oc.m ~ 40° where 
relativity reduces the nonrelativistic cross section by up to 5%. 
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FIG. 3: The relative deviation A from Eq. (|42|) for the three different relativistic approximations to the 
boosted potential V{fi)- For the description of the hnes see Fig|2] 



In ref. |1C| the nd total cross section has been investigated in a nonrelativistic scheme. Beyond 
that a very first step into relativity has been done using the optical theorem. In our notation it 
reads 



2 

atot = -^1'm ^ < (/) ;u„,/id|[/|(?:) ;U„,/id > . (43) 

Now we can investigate the changes in both ingredients on the right hand side due to relativ- 
ity, whereas in only the kinematical flux quantity |j| has been considered. The ratio of the 
relativistic to the nonrelativistic flux is given as 



where -E„ (-E'cz) is the neutron (deuteron) energy in the cm. system and gi^*^^'"^*^^ the relativistic or 
nonrelativistic relative momentum in the cm. system. In i.lO] only that ratio was considered which 
led to an increase in the total cross section by 3 (7) % at 100 (250 ) MeV. Now allowing also for a 
change of the nuclear matrix element < > (using the approximate boosted potential given 
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in Eq. (39)) the total cross section (not shown) is slightly smaller than the nonrelativistic one. In 
other words, the changes in < (f)\U\<j) > outweigh the kinematical effect, which by itself increases 
the total cross section. 

Now we come back to the elastic cross section which for the sake of completeness and clarity is 
shown. In our notation it has the form 

The kinematical factor in the bracket reduces to ^ in the nonrelativistic case. Again we can 
regard the ratio of the relativistic and nonrelativistic differential cross sections. They are of different 
type compared to the ratio of the total cross sections. Both ingredients, the kinematical factor 
and the nuclear matrix element enter now squared. It turned out as we have seen in Figs. 2 and 3 
that also in this case the dynamical effects caused by the decrease of the nuclear matrix element 
compensates the increase of the kinematical factor for most angles. Only at very backward angles 
a slight increase remains. 

Finally in Figsl^Hlwe compare relativistic and nonrelativistic predictions for the deuteron vector 
analyzing power iTn. Again the three approximations to the boosted potential V{q) are shown. It 
turns out that the relativistic effects are relatively small and they stay below ~ 5% in the angular 
regions outside of zero crossings. Other spin observables in elastic scattering behave similarly and 
are not shown. 



V. SUMMARY AND OUTLOOK 



We numerically solved the 3N Faddeev equation for nd scattering including relativistic features 
at the neutron lab energies E!^'^ = 28, 65, 135 and 250 MeV. The relativistic features are the 
relativistic form of the free propagator and the change of the NN potential caused by the boost of 
the 2N subsystem. In addition these boosts also induce Wigner spin rotations. For the momentum 
space basis we used the relative momentum of two free nucleons in their cm. system together with 
their total momentum which in the 3N cm. system is the negative momentum of the spectator 
nucleon. Such a choice of momenta is adequate for relativistic kinematics and allows to generalize 
the nonrelativistic approach used to solve the nonrelativistic 3N Faddeev equation to the relativis- 
tic case in a more or less straightforward manner. That relative momentum in the two-nucleon 
subsystem is a generalisation of the standard nonrelativistic Jacobi momentum p. The inclusion of 
the nucleon spins leads automatically to Wigner spin rotations in the context of boosting the 2N 
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FIG. 4: The deuteron vector analyzing power iTn in Nd elastic scattering. For description of lines see 
FigH The pd data at 28, 65, 135 and 250 MeV are from ref. S. and 

respectively. The pd 

data from l33| were taken at 190 MeV. 



cm. subsystems. The momentum partial wave basis, a generalisation of the nonrelativistic one, is, 
however, now more complex. As dynamical input we took the nonrelativistic NN potential AV18 
and generated in the 2N cm. system an exactly on-shell equivalent relativistic interaction v using 
an analytical scale transformation of momenta. We checked that in our energy range the boost 
effects for this potential could be sufficiently well incorporated by restricting the exact expression 
to the leading order terms in a q/uj and v/uj expansion. 

We found that in the studied energy range the effects of Wigner spin rotations are practically 
negligible for the cross section and analyzing powers. Relativistic effects for the cross section 
appear at higher energies and they are restricted only to the very backward angles where relativity 
increases the nonrelativistic cross section. At other angles the effects are small. In spite of the fact 
that the relativistic phase-space factor increases with energy faster than the nonrelativistic one, 
the relativistic nuclear matrix element outweighs this increase and leads for the cross section in 
a wide angular range to a relatively small relativistic effect. Also for spin observables (analyzing 
powers, spin correlation coefficients and spin transfer coefficients, not shown) no drastic changes 
due to relativity have been found. 

The comparison of our nonrelativistic theory with existing cross section data exhibits at the 
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FIG. 5: The same as in Figjjjbut for the deuteron vector analyzing power iTn. 

higher energies clear discrepancies. According to our presented results effects due to relativity 
are significant only in the region of very backward angles where they increase the cross section. 
They are relatively small in the region of the cross section minimum around 6 cm ~ 130°, where 
the discrepancies between the theory based on pairwise forces only and data are largest. At lower 
energies (up to about ~ 135 MeV) this discrep ancy can be removed when current three-nucleon 



forces (3NFs), mostly of 27r-exchange character 



33|, 3, 

are included in the nuclear Hamiltonian. 



At the higher energies, however, a significant part of the discrepancy remains and increases further 
with increasing energy. This indicates that additional 3N forces should be added to the 27r-exchange 
type forces. Natural candidates in the traditional meson-exchange picture are exchanges like n — p 



35l | in the order in which nonvanishing 3NF's 



and p — p- This has to be expected since in xPT 
appear the first time there are three topologies of forces, the 27r-exchange, a one-pion exchange 
between one nucleon and a two-nucleon contact interaction and a pure 3N contact interaction. 
They are of the same order and have to be kept together. Therefore it appears very worthwhile 
to persue a strategy adding in the traditional meson exchange picture further 3N forces. Our 
results here showing that relativistic effects based on relativistic kinematics and boost effects of 
the NN force are small support the usefulness of high energy elastic Nd scattering to study 3N 
force properties. 
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APPENDIX A: THE SPECIAL BOOST MATRIX (3{k) 

The special Lorentz transformation f3{k) = (3{ {ujm{k),k) ) is defined by f3{k){m,0) = k. It has 
the following matrix elements: 



Pi, = + — ' ' (i, j = 1, 2, 3) 

m[m + u}m\k)) 

= /3o^ = ^ (^ = 0,1,2,3). (Al) 

m 



APPENDIX B: THE MATRICES OF WIGNER ROTATION 

By straightforward calculation, the 3x3 matrix M related to the first equation in Eq. ()26() is 
presented. The four-momentum P is such that 

P{P) {Ujm{k),k) ) = Mp2),P2) 
f3{P) (LO^ik), -k) ) = {uJm{p3),P3)- (Bl) 

We obtain = 1,2,3) 

Mi,=6i, + fikikj + f2Pik, + hk,P, + hPiPj, (B2) 

where the four scalar functions /i, /2, /a and /4 depend on the momenta k and P: 

f ^ -Eq + Mq 

{m + ^) (k-P+^ + mMo) 



2 (4k- P + {Eo + Mo) {2m + Mq) 

f2 = ^ V 

(^0 + Mo) (2m + Mo) [2k-P + {Eq + 2m) Mo] 



-2 



2A:-P+(£;o + 2m) Mq 



(B5) 
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f 2m -Mo 

(^0 + Mo) (2 A; • P + (^0 + 2 m) Mo 



The quantities Mo and Eq are 



Mo = 2\]m? + k 2 (B7) 

and 



Eo = JM^ + P^. (B8) 



The matrix M related to the second equation in Eq. (|26|) results by replacing k by —A:. Clearly for 
P = there is no rotation. If one expands Mij in terms of the components Pi up to second order, 
a straightforward calculation leads to 



/3 



32 Sk^P^k-P , Mo-2m , ^^,-52, 



^ Mo2(2m + Mo)2 Mo2(2m + Mq) Mo2(2m + Mq) 

This is useful numerically in the determination of the Euler angles a, (3 and 7. For /? = only 
the combination (q + 7) occurs and one can put for instance a = 0. 

APPENDIX C: PERMUTATION OPERATOR 

Using Eq. (30) twice for the bra state 1 < kqa\ and the ket state \k'q'a' >2 one gets for the 
matrix element of the permutation operator in our partial wave basis: 



1 < k q a\ P \k' q a >i= 2 i <k q a \ k' q a >2= 

2 E E E E E 

mim2m3 ^l2^J.3^J.s ^J.l^l\^J.I^J. fi'^^i'^n^i ti^^i i-lj, i^l' 

(Xfix^rnillfij) {jnIfii\JM) {^f^2^f^3\sf^s) {lfiisfis\jf^) 

1 11 

iX'nx'-m2\I'ni') {j'fi'l'firlJM) {-fJ-2-f^3\s' fls') (^Vs'/is'liV) 

dq dq' 



N{q',-q-q') N{-q-q',q) 
S(k-\k{q',-q-q')\ ) 6{ k' - \k{ -q-q'^q] 



P k'^ 
Y.yiq) Yl;{k{ q',-q-q') ) Y^iq ') y/;,(fc( -q-q',q 

Z)i;,(P(/3(P( q',-q-q' )),k{ q',-q-q') )) 

Di*^,,{R{P{P{ q',-q-q' )),-k{ q'-q-q') )) 

Di,,>^{RW{P{ -q-q',q)),H -q-q',q))) 
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-q-q',q)),-k{ -q-q',q))) 
1 < ( (4)^ I >2, 



with 



H q',-q- q') = q" + ^q{l + yi{q,q' , x)) 
k{ -q-q',q) = - q- ^q"(.l + y2{q,q' ,x)). 

In Eq. (fn2|) occur 

yi{q,q',x) = 

E^, + Eg^+g^, + ^{E^, + E^+^,)^ - q^ 

with x = q - q', y2iq,q',x) = yi{q',q,x) and E^=u)m{q )• 

Expanding the product of the four D-matrices in Ea. ()Cl|l into spherical harmonics 

i?i;,(i?(/5(P( q',-q-q' )),k{ q',-q-q') )) 

^i^sWW q'-q-q' ))-k{ q',-q-q') )) 
Di^^,{RiP{P{ -q-q',q)),k{ -q-q\q))) 

-q-q',q)),-ki -q-q',q))) 

= E 4'I^l^'t:?'^''^'iq,q')yLAq)yL'M'iq'^ 

LML'M' 

and performing the integrations over q and q' in Ea. HCl|) leads to the following expression for 
matrix element of the permutation operator P: 

"1 , 5{k-ni) 5ik'-TT2) 



I < k q a\ P \k' q' a' >i = J dx 



^ -, — r Gaa'{q,q\x), 



Ni{q,q',x) N2{q,q',x) 
with 

G^Aq,q',x) = Y.Pk{x) E E "^^'^(1 + + 2/2)'^^'^'^'^'?, 

and 



vri = \jq''^ + ^q^{l + yif + qq'x{l+yi) 

71-2 = -sjq^ + \q''^0- + y2Y + qq'x{i + ^2) 



N,iq,q',x) ^ N{q'-q-q') 
N2iq,q',x) = N{-q-q',q). 
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The geometrical coefficients g'^^!^'^^^^{g,g') are given by 



with 



4vr \l/2 T t' 



(/')! 



{1)1 yTT' 



ci c' 



^{jfilfix + mi\Jn + fix + mi) (In- fj.2- M3S/^2 + fJ-sUfJ-) 

+ ^i\ + mi- fi'\jfi + fix + mi) (/'/i' - /"2 - ^35^2 + /^sIj'm') 

(AV + AtA + mi - ^' - m2-m2|/'/U + /xa + mi - /i') 

{L'M'X'fx + nx + mi - n' - m2\c[fj, + ^ua + mi - /x' - m2 + M') 

rr f2 fi I'] 

^(A;OriO|/iO)V/i(c'iO/iO|/iO) E(^0^20|/20) V/2(ci0/20|/20) ^ ^ 

^{limij2fJ' - fJ-2- ^J'3- mi^Hfj. - fj.2- fJ-s) 
mil 

(Ci - M - /Xa/2/U - /U2 - /U3 - + M + /iA|^2Ai - Ai2 - M3 - m^^ ) 

{fimfj2mf,yn' - ^3) 

(c'l^u + /iA + mi - /x' - m2 + M'fimf^\limi^). (C9) 

We used our standard notation 1 = 21 + 1. The coefficients ^^^^'^^^'"^'"^ {J'" aa' klil2l'il'2) are real 
and obey the following symmetry property 

^Z-MLri^"'^"™^"'"'('/'^««'fc/iZ2/'i/^) = A™i^?'"^"^^(J-aa'MiZ2Z'i/^). (CIO) 

The proof involves in addition to the definition ()C9|1 the properties /i + Z2 = ^, I'l +^2 = ^' aiid the 
parity conservation: (—1)'+'^ = (—1)''+'*''. The coefficients a'^jy/i/Mf"^^"^^ (q, q'), defined in HC4|) . 
are complex and obey 

*M2M3M2M3'^im2 / /N _ / sAf+M'+m2-At2-Ai3-M2+"^i-M'3+l -M2-/i3-/i2-M3-'»i-'»2 / /N /pii\ 
"'LML'M' K'djQ ) — \ ^) %_ML'-M' W^yj- V'-^J^J-j 

To obtain (ITTT]) one uses the relation Di^^,{a, f3,-f) = (-1)'^"'"' i^i^_„/(a, /5, 7)- The two last 
Clebsch-Gordan coefficients in Eq. HC9|) provide the condition that the phase factor in Eq. HC11|1 



is one. Thus the geometrical coefficients ^^'^/'^'^'^(g, (7') are real. 



23 



[1] R.B. Wiringa, V.G.J. Stoks, R. Schiavilla, Phys. Rev. C51, 38 (1995). 

[2] R. Machleidt, F. Sammarruca, and Y. Song, Phys. Rev. C53, R1483 (1996). 

[3] V.G.J. Stoks, R.A.M. Klomp, C.P.F. Terheggen, J.J. de Swart, Phys. Rev. C49, 2950 (1994). 

[4] J.L. Friar et al., Phys. Lett. B311, 4 (1993). 

[5] A. Nogga, D. Hiiber, H. Kamada, and W. Glockle, Phys. Lett. B409, 19 (1997). 

[6] H. Witala, W. Glockle, D. Hiiber, J. Golak, and H. Kamada, Phys. Rev. Lett. 81, 1183 (1998). 

[7] K. Sekiguchi et al., Phys. Rev. C65, 034003 (2002). 

[8] H. Witala et al., Phys. Rev. C63, 024007 (2001). 

[9] W.P. Abfalterer et al., Phys. Rev. Lett. 81, 57 (1998). 
[10] H. Witala et al., Phys. Rev. C59, 3035 (1999). 
[11] G. Rupp and J. A. Tjon, Phys. Rev. C45, 2133 (1992). 
[12] F. Sammarruca and R. Machleidt, Few-Body Syst. 24, 87 (1998). 

[13] A. Stadler, F. Gross, and M.Frank, Phys. Rev. C56, 2396 (1997); A. Stadler and F. Gross, Phys. Rev. 

Lett. 78, 26 (1997). 
[14] B. Bakamjian and L.H. Thomas, Phys. Rev. 92, 1300 (1952). 

[15] L.L. Foldy, Phys. Rev. 122, 275 (1961); R.A. Krajcik and L.L. Foldy, Phys. Rev. DIO, 1777 (1974). 
[16] B.D. Keister, in Few-Body Problems in Physics, edite by Paul Schoessow, AIP Conf. Proc. No. 334, 

(AIP, Woodbury, NY, 1995), p.l64. 
[17] B.D. Keister and W.N. Polyzou, Adv. Nucl. Phys. 20, 225 (1991). 
[18] H. Kamada, W. Glockle, J. Golak, and Ch. Elster, Phys. Rev. C66, 044010 (2002). 
[19] W. Glockle, T-S. H. Lee, and F. Coester, Phys. Rev. CSS, 709 (1986). 
[20] H. Kamada, W. Glockle, Phys. Rev. Lett. 80, 2547 (1998). 
[21] H.Witala, T.Cornelius and W.Glockle, Few-Body Syst. 3, 123 (1988). 
[22] W. Glockle, H.Witala, D.Hiiber, H.Kamada, J.Golak, Phys.Rep. 274, 107 (1996). 
[23] F. Coester, Helv. Phys. Acta 38, 7 (1965). 

[24] W.Glockle, The Quantum Mechanical Few-Body Problem, Springer- Verlag 1983. 

[25] S. Weinberg, The quantum theory of fields, vol. I, Cambridge University Press, 1995. 

[26] M.E. Rose, "Elementary theory of angular momentum" , Dover publications, Inc, New York, 1995. 

[27] H. Kamada, Few-Body Syst. Suppl.l2, 433 (2000). 

[28] K. Hatanaka et al., Nucl. Phys. A 426, 77 (1984). 

[29] S. Shimizu et al., Phys. Rev. C52, 1193 (1995). 

[30] K. Hatanaka et al., Phys. Rev. C66, 044002 (2002). 

[31] H.Witala et al., Few-Body Syst. 15, 67 (1993). 

[32] R.V. Cadman et al., Phys. Rev. Lett. 86, 967 (2001). 

[33] S.A. Coon et al., Nucl. Phys. A317, 242 (1979); S.A. Coon and W. Glockle, Phys. Rev. C23, 1790 



24 



(1981). 

[34] B.S. Pudliner et al., Phys. Rev. C56, 1720 (1997). 
[35] E. Epelbaum et al., Phys. Rev. C 66, 064001 (2002). 



25 



